Constraining Modified Gravity with Euclid 
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Future proposed satellite missions as Euclid can offer the opportunity to test general relativity on 
cosmic scales through mapping of the galaxy weak lensing signal. In this paper we forecast the ability 
of these experiments to constrain modified gravity scenarios as those predicted by scalar-tensor and 
f(R) theories. We found that Euclid will improve constraints expected from the PLANCK satellite 
on these modified gravity models by two orders of magnitude. We discuss parameter degeneracies 
and the possible biases introduced by modified gravity. 



I. INTRODUCTION 

Understanding the nature of the current observed ac- 
celerated expansion of our universe is probably the major 
goal of modern cosmology. Two possible mechanisms can 
be at work: either our Universe is described by general 
relativity (GR, hereafter) and its energy content is dom- 
inated by a negative pressure component, coined "dark 
energy" , either only " standard" forms of matter exist and 
GR is not valid on cosmic scales (see e.g. Q, @). 

All current cosmological data are consistent with the 
choice of a cosmological constant as dark energy compo- 
nent with equation of state w = P/p = — 1 where P and 
p are the dark energy pressure and density respectively 
(see e.g. i, i, 1). 

While deviations at the level of ~ 10% on w assumed 
as constant are still compatible with observations and 
bounds on w are even weaker if w is assumed to be 
redshift-dependent, it may well be that future measure- 
ments will be unable to significantly rule out the cosmo- 
logical constant value of w = — 1. 

Measuring w, however, is just part of the story. While 
the background expansion of the universe will be iden- 
tical to the one expected in the case of a cosmological 
constant, the growth of structures with time could be 
significantly different if GR is violated. Modified gravity 
models have recently been proposed where the expan- 
sion of the universe is identical to the one produced by 
a cosmological constant, but where the primordial per- 
turbations that will result in the large scale structures in 
the universe we observed today, grow at a different rate 
(see e.g. i, 0, 1). 

Weak lensing measurements offer the great opportu- 
nity to map the growth of perturbations since they re- 
late directly to the dark matter distribution and are not 
plagued by galaxy luminous bias (@, [l(|, [lljj). Re- 
cent works have indeed make use of current weak lensing 
measurements, combined with other cosmological observ- 
ables, to constrain modified gravity yielding no indica- 
tions for deviations from GR ([3, [ll, Qj, [H, QJ], 

E3). 

Next proposed satellite mission as Euclid ( 54], [l9j]) 
or WFIRST [20( could measure the galaxy weak lensing 



signal to high precision, providing a detailed history of 
structure formation and the possibility to test GR on 
cosmic scales. 

In this paper we study the ability of these future satel- 
lite missions to constrain modified gravity models and to 
possibly falsify a cosmological constant scenario. Respect 
to recent papers that have analyzed this possibility (e.g. 
[2l|, (22|) we improve on several aspects. First of all, we 
forecast the future constraints by making use of Monte 
Carlo simulations on synthetic realisations of datasets. 
Previous analyses (see e.g. Q, (26|, [27|) often used the 
Fisher matrix formalism that, while fast, it may lose its 
reliability when Gaussianity is not respected due, for in- 
stance, to strong parameter degeneracies. Secondly we 
properly include the future constraints achievable by the 
Planck satellite experiment, also considering CMB lens- 
ing, that is a sensitive probe of modified gravity (see e.g. 
[29f , [28[ and references therein) . Thirdly, we discuss the 
parameter degeneracies and the impact of modified grav- 
ity on the determination of cosmological parameters,. Fi- 
nally we focus on f(R) and scalar-tensor theories, using 
the general parametrization proposed by [111 . 

Our paper is structured as follows. In Section|TT]we in- 
troduce the parametrization used to describe departures 
from GR, and then specialize to the case of f(R) and 
scalar-tensor theories. In Section ILUl we describe Galaxy 
weak-lensing, while in section ITVl we discuss how to ex- 
tract lensing information from CMB data. We review the 
analysis method and the data forecasting in Section [V] 
In Section I VII we present our results and we derive our 
conclusions in Section IVIll 



II. MODIFIED GRAVITY PARAMETRIZATION 

In this section we describe the formalism we use to 
parametrize departures from general relativity. 

A. Background expansion 

In our analysis we fix the background expansion to a 
standard ACDM cosmological model. The reasons for 
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this choice are multiple; ACDM is currently the best fit 
to available data and popular models of modified grav- 
ity, e.g. f(R), closely mimic ACDM at the background 
level with differences which are typically smaller than the 
precision achievable with geometric tests [30(. The most 
significant departures happen at the level of growth of 
structure and, by restricting ourselves to ACDM back- 
grounds, we are able isolate them. 

B. Structure formation 



is modified with the introduction of a scalar field |32j. 
As shown in (2(| , for this kind of theories the parameters 
{/3j, Af } are related in the following way: 

ft = |=2-/3 2 | (5) 

and f < s < 4. 

This implies that we can analyze scalar-tensor theories 
adding 3 independent parameters to the standard cos- 
mological parameter set. 



In modified gravity models we expect departures from 
the standard growth of structure, even when the expan- 
sion history matches exactly the ACDM one. Dark mat- 
ter clustering, as well as the evolution of the metric po- 
tentials, is changed and can be scale-dependent. More- 
over, typically there might be an effective anisotropic 
stress introduced by the modifications and the two po- 
tentials appearing in the metric element, $ and ^, are 
not necessarily equal, as is in the ACDM model. Here 
we focus on the effect of the modified evolution of the 
potential, <f> + VP, on the CMB power spectra. 
In order to study the potentials evolution and 
to evaluate the growth of perturbations in mod- 
ified gravity models we employ the MGCAMB 
code developed in [26[ (and publicly available at 
http : //www . sf u . ca~g za5/MGCAMB . html| . In this 
code the modifications to the Poisson and anisotropy 
equations are parametrized by two functions p,(a, k) and 
7(0, k) defined by: 

a 2 

fc2 * = -^772-^(«,fc)pA , (I) 



2. f(R) theories 

In the specific case of f(R) theories we can additionally 
reduce the number of free parameters since f(R) theories 
correspond to a fixed coupling j3\ = 4/3 |33|. Moreover, 
to have ACDM background expansion the s parameter 
must be ~ 4 (2(|. The parametrization in Eq. (J3J effec- 
tively neglects a factor representing the rescaling of the 
Newton's constant (e.g. (l+fa) -1 in f(R) theories) that, 
as pointed out in [34| , is very close to unity in models that 
satisfy local tests of gravity [3(| and so negligible. How- 
ever, when studying the f(R) case, we need to include 
it to get a more precise MCMC analysis (see [34{ for the 
detailed expression of Eq. ([3])). Even with this extended 
parametrization, we have only one free parameter left, 
the length scale Ai. In this work we will constrain f(R) 
theories through this parameter, evaluating the effects of 
these theories on gravitational lensing. 



III. GALAXY WEAK LENSING 



"j- = 7(a, k) . 



(2) 



where pA = pS + 3^(p + P)v is the comoving density 
perturbation. These functions can be expressed u sing the 
parametrization introduced by [3fJ (and used in [26j): 



(i(a,k) = 



1 + k 2 a 
1 + Xj k 2 a s 



7(a,fc) = 



, + ftjAj k 2 a s 
1 + X\ k 2 a s 



(3) 



(4) 



where the parameters ft can be thought of as dimen- 
sionless couplings, Ai as dimensionful length scales and s 
is determined by the time evolution of the characteristic 
length scale of the theory. ACDM cosmology is recovered 
for /3i, 2 = f or X 2 2 = Mpc 2 . 



Being sensitive the the growth rate of the structure, 
weak lensing can be very useful to constrain modified 
gravity and to distinguish between various modified 
gravity models when combined with CMB observations. 
More generally, weak lensing (see 35] for a recent 
review or |http : / /w ww . gravitationallensing Tnetp 
is a particularly powerful probe for Cosmology, since 
it simultaneously measures the growth of structure 
through the matter power spectrum, and the geometry 
of the Universe through the lensing effect. Since weak 
lensing probes the dark matter power spectrum directly, 
it is not limited by any assumption about the galaxy 
bias (how galaxies are clustered with respect to the dark 
matter) that represents one of the main limitations of 
galaxy surveys. 

Following Q one can describe the distortion of the 
images of distant galaxies through the tensor: 



Scalar- Tens or theories tpij — 



-« - 7i -72 
-72 —K + 72 



This parametrization can be to constrain chameleon where k and (71,72) represents respectively the conver- 
type scalar-tensor theories, where the gravity Lagrangian gence (or magnification) and the shear (or stretching) 
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component of the distortion. The reconstruction of mat- 
ter density field can be conducted by looking at the corre- 
lations of the image distortions. The observable one has 
to deal with, will be hence a convergence power spectra 

(MHz!: 

P jk (t) = H$J o °°^W i (z)W j (z)P NL [P L (^y,*)] 

(6) 

where Pnl is the non-linear matter power spectrum at 
redshift z, obtained correcting the linear one Pl- W(z) 
is a weighting function, with subscripts i and j indicating 
the bins in redshifts. The function W(z) also encodes the 
cosmological information, being: 

Wi (z) = |o ro (i + Z ) T +1 (7) 

2 J Zz r(0,z') 

where: 

f dz> 

with E(z) = H(z)/ Hq and rii(z') is the fraction of sources 
belonging to the i — th bin. 

The observed convergence power spectra is affected 
mainly by a systematic arising from the intrinsic shear 
of galaxies rfms ■ This uncertainties can be reduced av- 
eraging over a large number of sources. The observed 
convergence power spectra will be hence: 

Cj k = Pjk + Sjkj^fij 1 (8) 

where fij is the number of sources per steradian in the 
j — th bin. 



IV. CMB LENSING EXTRACTION 

In the analysis we perform, we choose to introduce, in 
addition to galaxy weak lensing, the information derived 
from CMB lensing extraction. 

Gravitational CMB lensing, as already shown in Ref. [38[ , 
can improve significantly the CMB constraints on several 
cosmological parameters, since it is strongly connected 
with the growth of perturbations and gravitational po- 
tentials at redshifts z < 1 and, therefore, it can break 
important degeneracies. The lensing deflection field d 
can be related to the lensing potential 4> as d = V(j> [39| . 
In harmonic space, the deflection and lensing potential 
multipolcs follows: 

d? = -»v^+i)0r. (9) 

and therefore, the power spectra Cf d = (e?™d™*) and 
Cp = (<j>T<t>?*) are related through: 

Cf = i{l + \)Cf. (10) 



Gravitational lensing introduces a correlation between 
different CMB multipoles (that otherwise would be fully 
uncorrelated) through the relation: 

(afbf) = (-l) m S^5i'Cf + Y, Z?fL<t>L > (11) 

LM 

where a and b are the T, E, B modes and H is a linear 
combination of the unlensed power spectra C^ b (see 
for details). 

In order to obtain the deflection power spectrum from 
the observed C° b , we have to invert Eq. (fTTj) . defining a 
quadratic estimator for the deflection field given by: 

d(a,b)jf = nf W(a,b)??'f?a?bV' , (12) 

ii'rnm' 

where rif^ is a normalization factor needed to construct 
an unbiased estimator (d(a, b) must satisfy Eq. (|9])). This 
estimator has a variance: 

(d(o, &)f *d(a', &')£f > = 8? 5%' (Cf + Nf bb ') (13) 

that depends on the choice of the weighting factor W 
and leads to a noise iV£ a bb on the deflection power 
spectrum Cff obtained through this method. The choice 
of W and the particular lensing estimator we employ 
will be described in the next section. 



V. FUTURE DATA ANALYSIS 
A. Galaxy weak lensing data 



n ga i(arcmin 2 ) redshift f sky ^ ms 
35 < z < 2 0.5 0.22 



TABLE I. Specifications for the Euclid like survey consid- 
ered in this paper. The table shows the number of galaxies 
per square arcminute (n ga i), redshift range, f s t y and intrinsic 
shear (7™ s ). 



Future weak lensing surveys will measure photometric 
redshifts of billions of galaxies allowing the possibility of 
3D weak lensing analysis (e.g. |4lU44j )) or a tomographic 
reconstruction of growth of structures as a function of 
time through a binning of the redshift distribution of 
galaxies, with a considerable gain of cosmological infor- 
mation (e.g. on neutrinos [Ha]; dark energy the 
growth of structure [3 [4?| and map the dark matter 
distribution as a function of redshift (48J). 
Here we use typical specifications for futures weak lens- 
ing surveys like the Euclid experiment, observing about 
35 galaxies per square arcminute in the redshift range 
< z < 2 with an uncertainty of about a z = 0.03(1 + z) 
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(see [13), to build a mock dataset of convergence power 
spectra. Table U shows the number of galaxies per 
arcminute -2 (n ga i), redshift range, f s k v and intrinsic 
shear for this survey. The expected la uncertainty on 
the convergence power spectra P(£) is given by (49|: 



{21 + l)f sky A e 



P{£) 



7 2 

irms 
n g al 



(14) 



For the convergence power spectra we use l ma x = 1500 in 
order to exclude the scales where the non-linear growth 
of structure is more relevant and the shape of the non- 
linear matter power spectra is, as a consequence, more 
uncertain (see [52]). We calculate the power spectra at 
a mean redshift z — 1. Ag in the (fT4|) is the bin used 
to generate data. Here we choose Ai = 1 for the range 
2 < I < 100 and A e = 40 for 100 < I < 1500. 
In this first-order analysis we are not considering other 
systematic effects as intrinsic alignments of galaxies, se- 
lection effects and shear measurements errors due to un- 
certainties in the point spread function (PSF) determi- 
nation. Of course future real data analysis will require 
the complete treatment of these effects in order to avoid 
biases on the cosmological parameters. 



B. CMB data 

We create a full mock CMB datasets (temperature, 
E-polarization mode and lcnsing deflection field) with 
noise properties consistent with the Planck 50] experi- 
ment (see Tab. HTlfor specifications). 

Experiment Channel FWHM AT /T 



Planck 70 


14' 


4.7 


100 


10' 


2.5 


143 


7.1' 


2.2 


fsky = 0.85 







TABLE II. Planck experimental specifications. Channel 
frequency is given in GHz, FWHM (Full- Width at Half- 
Maximum) in arc-minutes, and the temperature sensitivity 
per pixel in jj,K/K. The polarization sensitivity is AE/E = 
AB/B = V2AT/T. 



We consider for each channel a detector noise of 
ur 1 = (9a) 2 , where 9 is the FWHM (Full- Width at Half- 
Maximum) of the beam assuming a Gaussian profile and 
a is the temperature sensitivity AT (see Tab. UU for the 
polarization sensitivity). We therefore add to each Cg 
fiducial spectra a noise spectrum given by: 



N i = w- 1 exp(£(£ + l)/£l), 
where £{, is given by it, = \/8\n2/9. 



(15) 



In this work, we use the method presented in [40j to 
construct the weighting factor W of Eq. (TT21 . In that pa- 
per, the authors choose IF to be a function of the power 



spectra C^ b , which include both CMB lensing and pri- 
mary anisotropy contributions. This choice leads to five 
quadratic estimators, with ab = TT,TE,EE,EB,TB; 
the BB case is excluded because the method of Ref. 
is only valid when the lensing contribution is negligible 
compared to the primary anisotropy, assumption that 
fails for the B modes in the case of Planck. 
The five quadratic estimators can be combined into a 
minimum variance estimator which provides the noise on 
the deflection field power spectrum Cf d : 



n; 



dd 



1 



aba'b' \ — l 



(16) 



We compute the 
ing noise for 



minimum variance lens- 
Planck experiment by 
means of a routine publicly available at 
http : //lesgourg . web . cern . ch/lesgourg/ codes . html. 
The datasets (which include the lensing deflection power 
spectrum) are analyzed with a full-sky exact likelihood 
routine available at the same URL. 



C. Analysis method 

In this paper we perform two different analysis. First, 
we evaluate the achievable constraints on the f{R) 
parameter A 2 and on the more general scalar-tensor 
parametrization including also /?i and s. Secondly, we 
investigate the effects of a wrong assumption about the 
modified gravity on the cosmological parameters, by gen- 
erating an f(R) datasets with non-zero A 2 fiducial value 
and analysing it fixing \\ — Mpc 2 . We conduct a full 
Monte Carlo Markov Chain analysis based on the pub- 
licly available package cosmomc [53] with a convergence 
diagnostic using the Gelman and Rubin statistics. 
We sample the following set of cosmological parameters, 
adopting flat priors on them: the baryon and cold dark 
matter densities flbh 2 and £l c h 2 , the ratio of the sound 
horizon to the angular diameter distance at decoupling 
9 S , the scalar spectral index n s , the overall normaliza- 
tion of the spectrum A s at k = 0.002 Mpc -1 , the optical 
depth to reionization t, and, finally, the modified gravity 
parameters X\ , /3i and s. 

The fiducial model for the standard cosmological param- 
eters is the best-fit from the WMAP seven years anal- 
ysis of Ref. [H| with n b h 2 = 0.02258, fl c h 2 = 0.1109, 
n s = 0.963, t = 0.088, A s = 2.43 x 10 -9 , 6 = 1.0388. 
For modified gravity parameters, we first assume a fidu- 
cial value \\ = Mpc 2 and fix /3i = 1.33 and s = 4 
to test the constraints achievable on the f{R) model. 
We then repeat the analysis allowing /3i and s to vary. 
Furthermore, to investigate the ability of the combina- 
tion of Planck and Euclid data to detect an hypothet- 
ical modified gravity scenario, we study a model with 
fiducial \\ = 300 Mpc 2 leaving A 2 , f3\ and s as free 
variable parameters allowing them to vary in the ranges 
< X\ < 10 6 , 0.1 < j3x < 2 and 1 < s < 4. Finally, we 
analyse a dataset with a fiducial value A 2 = 300 Mpc 2 
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but assuming a wrong ACDM model fixing A 2 = Mpc 2 , 
to investigate the bias introduced on the cosmological 
parameter due to a wrong assumption about the gravity 
model. 



VI. RESULTS 

In Table Hill we show the MCMC constraints at 68% c.l. 
for the f(R) case for Planck alone and Planck combined 
with Euclid. For this last case we also fit the data fixing 
A 2 to 0, thus performing a standard analysis in a General 
Relativity framework, in order to show the importance of 
the degeneracies introduced by \\ on the other cosmolog- 
ical parameters errors. The parameters mostly correlated 
with modified gravity are Hq and Q, c h? (see also Figure[T]) 
because these parameters strongly affect the lensing con- 
vergence power spectrum as well as \\ through P(k,z). 
As expected in fact, when assuming general relativity we 
find strong improvements on the errors on these parame- 
ters for the combination Planck+Euclid in comparison to 
the varying A 2 analysis. We note that the constraints on 
the standard cosmological parameters are in good agree- 
ment with those showed in 15411. 



effect of the modified gravity on the growth of structures. 





Planck 


Planck+Euclid 


Fiducial: 

Model: 

Parameter 


xi = 

varying A x 


A? = 
varying A? 


X{ = o 

fixed A? 


A(fW) 

A(n c h 2 ) 

A(9 S ) 

A(r) 

A(n s ) 

A(log[10 10 A s ]) 
A(H ) 
A(Ov) 
A 2 (Mpc 2 ) 


0.00013 
0.0010 
0.00027 
0.0041 
0.0031 
0.013 
0.50 
0.0050 
< 2.42 x 10 4 


0.00011 
0.00073 
0.00025 
0.0030 
0.0029 
0.0091 
0.38 
0.0040 
< 2.9 x 10 2 


0.00010 
0.00057 
0.00023 
0.0026 
0.0027 
0.0091 
0.29 
0.0031 



TABLE III. 68% c.l. errors on cosmological parameters. Up- 
per limits on A 2 are 95% c.l. constraints. 



In the third column we show constraints on the cosmological 
parameters when fitting the data assuming general relativity, 
i.e. fixing A 2 = Mpc 2 . 



In Figure [T] we show the 68% and 95% confidence 
level 2-D likelihood contour plots in the ft m — A?, 
Hq — \\ and n s — \\ planes, for Planck on the left 
(blue) and Planck+Euclid on the right (red). As 
one can see the inclusion of Euclid data can improve 
constraints on the standard cosmological parameters 
from a 10% to a 30%, with the most important im- 
provements on the dark matter physical density and 
the Hubble parameter to which the weak lensing is of 
course very sensitive as showed by Eq. (|6|) and (J7J. 
Concerning modified gravity, Euclid data are decisive 
to constrain A 2 , improving of two order of magnitude 
the 95% c.l. upper limit, thanks to the characteristic 





Planck+Euclid 


Planck+Euclid 


Fiducial values 


Model: 


Xi = o 


varying A^ 




Parameter 










0.022326 ± 0.000096 


0.02259 ± 0.00012 


0.02258 


n c h 2 


0.1126 + 0.00055 


0.11030 ±0.00083 


0.1109 


0s 


1.0392 + 0.00023 


1.0395 ± 0.00025 


1.0396 


T 


0.0775 ± 0.0024 


0.08731 ± 0.0029 


0.088 


n s 


0.9592 ± 0.0027 


0.9636 ± 0.0029 


0.963 


Ho 


69.94 + 0.27 


71.20 + 0.42 


71.0 


n A 


0.724 ± 0.003 


0.738 ± 0.005 


0.735 


C"8 


0.8034 ± 0.0008 


0.8245 ± 0.0039 


0.8239 



TABLE IV. best fit value and 68% c.l. errors on cosmological 
parameters for the case with a fiducial model A 2 = 300 fitted 
with a ACDM model where A 2 = is assumed. 



Moreover, when analyzing the f(R) mock datasets 
with \\ — 300 Mpc 2 as fiducial model, assuming A 2 — 
Mpc 2 we found a consistent bias in the recovered best 
fit value of the cosmological parameters due to the 
degeneracies between \\ and the other parameters. 
As it can be seen from the comparison of Figures 
Q] and Figures [2] and from table IIVI the shift in the 
best fit values is, as expected, along the degeneracy 
direction of the parameters with A 2 , for example for 
n s , Hq and £l m - These results show that for an even 
small modified gravity, the best fit values recovered 
by wrongly assuming general relativity are more than 
68% c.l. (for some parameters at more than 95% c.l.) 
away from the correct fiducial values, and may cause an 
underestimation of n s and Hq and an overestimation of 
<7g and f2 m . More generally, as shown in table IIVI all 
parameters are affected. 

We conclude, hence, that a future analysis of so high 
precision data from Euclid and Planck will necessarily 
require to allow for possible deviations from general 
relativity, in order to not bias the best fit value of the 
cosmological parameters. 

We also perform an analysis allowing ft and s to vary; 
in this way we can constrain not only /(i?)theories but 
also more general scalar-tensor models, adding to the 
standard parameter set the time variation of the new 
gravitational interaction s and the coupling with matter 
ft- 

We perform this analysis assuming as a fiducial model a 
f(R) theory with A 2 = 3.0 x 10 4 Mpc 2 and ft = 4/3. 

In Table [V] we report the 68% c.l. errors on the 
standard cosmological parameters, plus the coupling 
parameter ft. Performing a linear analysis, with a 
fiducial value of A 2 = 3 x 10 4 , we obtain constraints 
on ft with A(ft) = 0.038 at 68% c.l. and therefore 
potentially discriminating between modified gravity 
models and excluding the ft — 1 case (corresponding to 
the standard ACDM model) at more than 5 — a from 
a combination of Planck+Euclid data (only 2 — a for 
Planck alone). 
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Planck 


Planck+Euclid 


Fiducial: 


A? = 3.0 x 10 4 


A? = 3.0 x 10 4 


Parameter 






A(Q b h A ) 


0.00013 


0.00011 


A(Q. c h 2 ) 


0.0011 


0.00082 


A(9 a ) 


0.00026 


0.00025 


A(r) 


0.0043 


0.0040 


A(n s ) 


0.0033 


0.0029 


A(log[10 10 A s ]) 


0.014 


0.011 


A(ifo) 


0.54 


0.40 


A(fi A ) 


0.0060 


0.0045 


A(/3i) 
A? 


0.13 


0.038 


unconstrained 


unconstrained 


s 


unconstrained 


unconstrained 



TABLE V. 68% c.l. errors on cosmological parameters and /3i. We do not show limits on Ai and s because this kind of analysis 
does not allow to constrain them (see text). 



The strong correlation present between /3i and \\ (see 
cq. [3]) implies that, choosing a lower \\ fiducial value for 
a f(R) model, the same variation of /?i brings to smaller 
modifications of CMB power spectra and therefore we 
can expect weaker bounds on the coupling parameter. 
In order to verify this behaviour we made three analysis 
fixing s = 4 and choosing three different fiducial values 
for A?: 3 x 10 2 , 3 x 10 3 and 3 x 10 4 Mpc 2 . The 
respectively obtained Pi 68% c.l. errors are 0.11, 0.052 
and 0.035, confirming the decreasing expected accuracy 
on j3i for smaller fiducial values of A 2 . 

The future constraints presented in this paper are ob- 
tained using a MCMC approach. Since most of the fore- 
casts present in literature on f(R) theories are obtained 
using a Fisher matrix analysis, it is useful to compare 
our results with those predicted by a Fisher Matrix ap- 
proach. We therefore perform a Fisher Matrix analysis 
for Planck and Planck+Euclid (see [55T - I571 ] ) assuming a 
ACDM fiducial model and we compare the results with 
those in Table Gill 

We find that for Planck alone the error on Ai is under- 
estimated by a factor ~ 3 while the error is closer to 
the MCMC result for the Planck+Euclid case (underes- 
timated by a factor ~ 1.2). 

VII. CONCLUSIONS 

In this paper we forecasted the ability of future weak 
lensing surveys as Euclid to constrain modified gravity. 



We restricted our analysis to models that could mimic a 
cosmological constant in the expansion of the Universe 
and can therefore be discriminated by only looking at 
the growth of perturbations. We have found that Eu- 
clid could improve the constraints on these models by 
nearly two order of magnitudes respect to the constraints 
achievable by the Planck CMB satellite alone. We have 
also discussed the degeneracies among the parameters 
and we found that neglecting the possibility of modified 
gravity can strongly affect the constraints from Euclid on 
parameters as the Hubble constant Hq, fi m and the am- 
plitude of r.m.s. fluctuations erg. In this paper we found 
that, considering more general expansion histories, would 
further relax our constraints and increase the degenera- 
cies between the parameters. However other observables 
can be considered as Baryonic Acoustic Oscillation and 
luminosity distances of high redshift supernovae to fur- 
ther probe the value of w and its redshift dependence. 
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FIG. 1. 2-dimensional contour plots showing the degeneracies at 68% and 95% confidence levels for Planck on the left (blue 
countours) and Planck+Euclid on the right (red countours). Notice different scale for abscissae. 




FIG. 2. 2-dimensional contour plots showing the degeneracies at 68% and 95% confidence levels for Planck+Euclid assuming a 
f(R) fiducial cosmology with Ai = 300Mpc 2 considering an analysis with Ai fixed to (blue countours) or allowing it to vary 
(red countours). 



